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Abstract 

The unfolding of a gamma ray spectrum experience many difficulties due to 
noise in the recorded data, that is based mainly on the change of photon 
energy due to scattering mechanisms (either in the detector or the medium) , 
the accumulation of recorded counts in a fixed energy interval (the channel 
width of the detector) and finally the statistical fluctuation inside the detec- 
tor. In order to deal with these problems, a new method is developed which 
interpolates the ideal spectrum with the use of special designed derivative 
kernels. Preliminary simulation results are presented and show that this 
approach is very effective even in spectra with low statistics. 
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1. Introduction 

Response function calculation of Nal based scintillators has many appli- 
cations like process control tasks in manufacturing industry, in oil detection, 
in safety and alarm systems, in Prompt Gamma Neutron Activation Analy- 
sis (PGNAA) and others (see for example [H,0], 0] and references there in). 
In general, Monte Carlo techniques are used to calculate the interactions of 
source photons with the detector (Q, j^], 0]) and thus the response function. 
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Also both analytical ([7|, [8() and statistical ([9( techniques have been used 
too. 

On the other hand, during the detection of gamma rays, several problems 
are encountered, ie. the efficiency vs. resolution of semiconductor or scin- 
tillation detectors used, the geometry, which causes in turn uncertainties in 
the solid-angle determination, etc. Also the form of the spectrum becomes 
more complex due to the following properties: (i) the scintillation detectors 
have a lower energy resolution compared to Ge detectors, (ii) the environ- 
ment and/or shielding play an important role because of the scattering of 
high energy x-rays into the detector, (iii) the first and second escape peaks 
become important at high energies and (iv) a significant tail develops to- 
wards the low-energy continuum due to Compton scattering and escape of 
bremsstrahlung from the detector. These effects reduce the detection effi- 
ciency in the full-energy peak, and have also other serious consequences. If 
the spectrum is complex, with a continuous 7-yield (e.g. due to statistical 
7-rays following the decay of highly excited nuclei), the large superimposing 
continuous tails of the high-energy 7-rays may hamper an accurate evaluation 
of the continuous 7-yield. 

To improve these drawbacks several attempts have been made in the past. 
In the experiments a combination of different detectors (Ge and BaF, anti- 
Compton shields, etc.) has been used. However, these techniques either 
reduce the overall efficiency by rejecting a large part of the detected events 
(anti- Compton), or hamper a precise determination of the overall efficiency 
(addition of coincident signals from different types of detectors). In the data 
analysis the generally used forward method fits the measured spectrum using 
appropriate physical models (input information): a master-spectrum is gen- 
erated using e.g. statistical model calculations (some model parameters are 
to be adjusted later), which is then folded with the detector response function 
and the resulting spectrum is compared with observation. Finally, the model 
parameters are adjusted, until an acceptable agreement is found. Problems 
arise here from peaks in the experimental spectrum due to contaminants in 
the target creating discrete lines, which cannot be simulated easily. The 
remaining problem is the choice of the physical model and the appropriate 
model parameters. If several physical processes compete, the generation of 
the master spectrum can often be ambiguous (Io| . On the other hand, even if 
the model spectrum is accurate, the accuracy of the unfolding process is re- 
duced due to two main reasons: (i) the noise in the measuring spectrum and 
(ii) the fact that the measuring spectrum represents the total counts recorder 
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in a finite energy interval, which is the channel width of the detector. 

The purpose of this work is to present a new method which can improve 
the unfolding procedure of a given measured spectrum. The method interpo- 
lates the ideal spectrum with the use of special designed derivative kernels. 
Preliminary simulation results are presented which show that this approach 
is very effective even in spectra with low statistics. 

2. Derivative kernels in unfolding procedure 

Consider the case where a radioactive source emits photons in a uniform 
medium and at a given point a photon detector has been placed. Photons, 
after their emission and before they reach the detector, interact with the 
atoms of the uniform medium and can change their energy due to Compton 
scattering or pair production, or disappear due to the photoelectric effect. 
The effect of the interaction of photons with the medium can be formulated 
as follows. Let S(E) be the source spectrum and M(E) the measured one. 
In vacuum, 



where R(E, V) is equal to the number of photons that will be recorded at 
energy E when one photon is emitted with energy V. The function R(E, V) 
is known as the transfer function of the detector. In the uniform medium, this 
relation is more complicated. If a photon with initial energy U is emitted, 
then there is a probability P(V, U) that the photon will reach the detector 
surface with a final energy V. Thus, the measured spectrum now will be 
given by: 



can be regarded now as the modified transfer function of the detector, for 
operation inside the uniform medium [11]. The measured spectrum can now 
be expressed as: 




(1) 





Changing the order of integration, the function 




(3) 





3 



But instead of the function M(E) the detector integrates this function in 
small energy intervals, called channels. Thus, the detector output M(E) is 
given: 

pE+e 

M(E) = / M(V) ■ dV (5) 
Je 

where e is the channel width. Consider now the function 

pE+E' 

m(E,E') = / M(V) ■ dV (6) 
Je 

Since M(E) is equal to m(E,e), the 2-dimensional function m(E,E') is 
known on the grid (n\ ■ e, ni • e), n\,n2 = 0(1) N. Our purpose now is to 
find optimal derivative kernels in order to calculate derivatives of the func- 
tion m(E,E'). Then, we can calculate M(E): 

= M(E) (7) 

An important property of m(E, E') which allows for the application of equa- 
tion is that: 

m(E,-E') = -m(E-E',E') (8) 
Figure [T] shows the measured spectrum from a Nal detector in an underwater 



experiment, described in [12|]. Both M(E) and the calculated from equation 
([70 M(E) are shown in Figure [2J 

The unfolding of the gamma ray spectrum M(E) can now be easily ob- 
tained in the case where radioactive sources emit photons in discrete energies 
and the counting rate is low enough to avoid additive effects in the detector. 
In this case and based on the linearity of the folding mechanism, we assume 
that 

k 

S(E) = Y,an5(E-E n ) (9) 

n=0 

and we want to calculate both a n and E n . Then, it is easily found that 

k 

M(E) = J2^R(E,E n ) (10) 

n=l 

Finally, consider a continuous function g : R 2 — » R and its discrete version 

oo 

9s = ^2 g(x,y)5(x-n l T)-5(y-n 2 T) (11) 

ni,n2=—oo 
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where 5 is the Dirac delta function. The knowledge of the discrete version g s 
can lead to the reconstruction of the continuous function g with the aid of a 
kernel K such that 

oo 

f}fay) = ^2 9s(ri 1 ,n 2 )K(x-n 1 T,y-n 2 T) (12) 

711,712=— OO 

The ideal interpolation where g — g is achieved if K(x,y) = st(x) ■ sr(y), 
where 

= ^ 

7ra;/i 

and T is the Nyquist rate. For practical reason, we assume that 

K(x,y) = d (x)-d Q (y) (14) 

and 

^ = *.(*) («» 

Then, the expression for the derivative with respect to x of the reconstructed 
function g becomes: 

oo 

D x {g}(x,y)= ^ gs(n 1 ,n 2 )d 1 (x-n 1 T)d (y-n 2 T) (16) 

m ,7i2=— oo 

In order to construct an efficient kernel, it is not necessary that di(x) = d' (x). 
Although this seems controversial consider the following example: it is com- 
mon to use a sampled Gaussian and its derivative. However, because the 
Gaussian is not strictly bandlimited, sampling introduces artifacts, thus de- 
stroying the derivative relationship between the resulting kernels. So, instead 
we choose to simultaneously design a pair of discrete kernels that optimally 
preserve the required derivative relationship. If 

D (co) = J2do(nT)e-^ n , AH^^Kje"™ (17) 

n n 

with u = 2n /T are the discrete Fourier transforms of do, d\, we can demand 
that 

icuDoito) = Di(w) (18) 

in the case of one dimensional signals g. In the case of two dimensional 
signals, we can demand for example that the pair of kernels preserve the 



derivative relationship in all directions 13 
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3. Simulation results 



In order to test the new method, a simulation experiment was performed. 
A folded spectrum is produced using the transfer function for a Nal based 



measuring system calculated in [141 ] . The spectrum is folded again using 
the method presented in 15( in order to simulate an underwater measuring 
system. Several simulated spectra were produced, with different number of 
photo-peaks and different number of total recorded counts in order to account 
for the spectrum statistics. Figure [3] shows the overall error in the unfolded 
spectrum, using a Gaussian derivative kernel (o) and three derivative kernels 
DK3 (□), DK4 (x) and DK5 (o) with 3,4 and 5 points respectively calculated 
in [l3| . Furthermore, in Figure H] the dependence of the error in the unfolded 
spectrum on spectrum statistics is shown for the Gaussian (o) and DK5 (x) 
derivative kernel. It is clear that the new method is very promising even in 
cases with low statistics. A special experimental setup is under construction 
to test the new method in real spectrums. Moreover, new derivative kernels 
are designed in order to optimize their behavior. 



4. Conclusions 

Preliminary results on interpolation of a measured spectrum with deriva- 
tive kernels, show that the unfolding procedure becomes more accurate even 
in cases of low statistics. The use of derivative kernels facilitate the numer- 
ical differentiation which is of high importance in both peak detection and 
spectrum unfolding. 
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Figure 1: Two dimensional representation of measured spectrum obtained from [l 
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Figure 3: Unfolding of simulated spectra generated with the method in DK3 uses the 
3-point kernel couple, DK4 the 4-point one and DK5 the 5-point kernel. The Gaussian 
curve is produced using a sampled Gaussian and its derivative. Gaussian derivative kernel 
(o), DK3 (□), DK4 (x) and DK5 (o). 
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Figure 4: Dependence of unfolding accuracy on spectrum statistics for the Gaussian (o) 
and DK5 (+) case. The horizontal axis represents the total number of counts recorded. 
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